Compact stars in Eddington inspired gravity 
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A new, Eddington inspired theory of gravity was recently proposed by Banados and Ferreira. 
It is equivalent to General Relativity in vacuum, but differs from it inside matter. This viable, 
one parameter theory was shown to avoid cosmological singularities and turns out to lead to many 
other exciting new features that we report here. First, for a positive coupling parameter, the field 
equations have a dramatic impact on the collapse of dust, and do not lead to singularities. We 
further find that the theory supports stable, compact pressureless stars made of perfect fluid, which 
provide interesting models of self-gravitating dark matter. Finally, we show that the mere existence 
of relativistic stars imposes a strong, near optimal constraint on the coupling parameter, which can 
even be improved by observations of the moment of inertia of the double pulsar. 
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I. Introduction. Einstein's General Relativity (GR) 
is able to explain a wide variety of phenomena at solar 
system scales and beyond, and after decades of intense 
scrutiny stands as the most attractive theory of gravity. 
However, high-curvature corrections may be necessary to 
address unresolved issues, such as the presence of singu- 
larities in cosmology and in the interior of black holes. 

Recently, an intriguing alternative to GR (based on 
an original proposal by Eddington) was put forward by 
Banados and Ferreira (BF) [l[ (see also @). BF theory 
is described by the action 



S=^-Jd 4 x (V- \9ab + KR ab (T)\ - Ay^s) , 



(1) 



where R ab (T) denotes the symmetric part of the Ricci 
tensor, built from the connection T c ab and A is related to 
the cosmological constant, A = (A — 1)/k. We will focus 
on asymptotically flat solutions and set A = 1. The met- 
ric g and the connection T are independent fields and, 
at the classical level, matter is minimally coupled to the 
metric only. The BF proposal explores the fact that the 
coupling between matter and gravity is one of the least 
tested sectors of gravitation; in fact the theory can be 
shown to be completely equivalent to GR in vacuum [lj. 
However, it dramatically differs from GR in the pres- 
ence of matter, for example it yields a singularity-free 
cosmology, thus presenting itself as a potentially excit- 
ing gravity theory. The BF theory modifies the Newto- 
nian regime (again, in the matter coupling), but tests of 
gravity within matter are extremely hard to carry, partly 
because we understand the coupling to matter so poorly. 
Thus, constraints on the BF theory as a result of Earth- 
based experiments are hard to accomplish. 

Here we show that non-linear effects and deviations 
from Einstein's theory are more pronounced inside high- 
density objects. When nR a b <C 1, the first corrections to 
the Einstein equations read 



Rab(T) = T ab - -jTg ab + K 



Sab 



■^Sg a b 



where two types of 0{n) corrections appear: those hidden 
in i? a f,(r), which implicitly depend on derivatives of mat- 
ter fields, and those depending on S a b = T c a T c b~ ^TT a b, 
which arc quadratic in the matter fields. Hence we ex- 
pect strong corrections at high densities or where strong 
matter gradients exist, for example in early cosmology [l| 
or inside neutron stars (NSs). Higher order corrections 
in the matter fields were also discussed in Ref. [1] to cure 
cosmological singularities. The purpose of this letter is 
to show that the best possible constraints on the theory 
arise from the study of NSs and other compact objects. 
In the process, we report new remarkable features of BF 
theory. 

II A . Stars in the non-relativistic limit. Let us start 
by discussing the nonrelativistic limit of ([TJ. The modi- 
fied Poisson equation reads [l[ 



V 2 $ = 4irGp + nV 2 p/4. 



(2) 



From Eq. ([2]) , and requiring spherical symmetry, the hy- 
drostatic equilibrium equation follows 



dP/dr = -Gm{r)p/r 2 - Kpp'/A. 



(3) 



While BF corrections are absent for constant density pro- 
files, interesting effects may show up for non-trivial mat- 
ter distributions. Newtonian stellar models are solutions 
of Eq. supplemented by the standard mass conserva- 
tion, dm/dr = 47tr 2 p(r) and an equation of state (EOS). 
We note that constant density stars in BF theory are po- 
tentially pathological, since they introduce a Dirac delta 
contribution in Eq. ©. For this reason, in this letter 
we shall focus on more realistic, polytropic models of the 
form P(p) = Kp( n+1 " n , where K and n are constants. 
IIAl. Newtonian pressureless stars. Remarkably, this 
theory supports pressureless stars, i.e. stars made of non- 
interacting particles, which provide interesting models for 
self-gravitating dark matter. Indeed, if P = and k > 0, 
Eq. © is solved by 



+ 0(k 2 ), 



p(r) = p c sm(wr) / (tur) , 



VGtt/k. (4) 



2 



The radius and mass of the star read R = -k/td and 
M = 47r 2 p c /-n7 3 , respectively. In the interior, the New- 
tonian potential is constant and it matches continuously 
the vacuum potential M / r at the radius. Below, we prove 
that these solutions are also linearly stable. 
IIA2. Newtonian polytropic models. For a generic poly- 
tropic index n, the held equation must be solved numer- 
ically, imposing p ~ p c + p 2 r 2 at the center. It is easy 
to show that realistic stellar configurations (with p->0 
at the surface of the star) can only exist provided the 
following condition is satisfied 



k > 



-4K(1 + l/n)p, 



-l+i/i 



(5) 




t/ff 



FIG. 1. Oscillatory behavior of collapse-related quantities as 
functions of time for k\t)\ = 6.5. 



Similar constraints exist for any EOS for which the pres- 
sure increases monotonically with the density. For k > 0, 
condition ((SJ) is always fulfilled. In some cases the Lane- 
Emdcn equation obtained from Eq. ([3]) can be solved an- 
alytically Q. For instance if n = I, P(p) = Kp 2 , and the 
solution reads as in 0, but with ru = i^Gn/ (8K + k), 
so that it exists for n > —8K and reduces to the pres- 
surclcss case for K = 0. 

IIB. Stability in the non-relativistic limit. We now 

discuss stability of the Newtonian configurations against 
radial perturbations. The standard treatment can be ex- 
tended straightforwardly to encompass BF theory Q. 
Assuming a time dependence ~ e lut for the fields, the 
modified eigenvalue equation reads 



4£P' up 
r 4 



-tp'-i'p' 
r 



where 7 is the adiabatic index of the perturbations. This 
equation must be solved for the Lagrangian displacement 
£ requiring regularity at the center and at the radius. An 
instability corresponds to an eigenmode with cj 2 < 0. 
IIB1. Pressureless stars. For P = and p given in 
our numerical study found eigenmodes with uj 2 > and 
no unstable mode. Therefore, pressureless stars in the 
modified Newtonian theory are stable. As we shall see, 
these solutions persist in the fully relativistic theory. 
IIB2. Newtonian polytropic stars. In Newtonian gravity, 
polytropic models with 7 = 4/3 are marginally stable 
for any polytropic index n [f|. In our case, these models 
are stable if K > and unstable if n < 0. For generic 
values of 7, positive values of k contribute to stabilize 
the models, while negative values work in the opposite 
direction. 

IIC. Gravitational collapse. The collapse of incoher- 
ent dust in the Newtonian limit shares many properties 
with its relativistic analogue 0, ■ The relevant Eulerian 
equations governing the fluid dynamics are Eqs. (9) and 
(f I) in Ref. together with 

cVt(x) + u(x)0 r it(x) = -GM(x)/r 2 - Kd r p(x)/4 , 

where w(x) is the fluid velocity and x = (t, r). In stan- 
dard Newtonian gravity, k = 0, the equations can be 
solved analytically when /o(x) = pit) and they correspond 



to the relativistic Oppcnheimer-Snyder solution Q. The 
dust collapses in a finite time tc = tt y/ R 3 / \8Mt) , where 
R and Mr are the initial radius and the total mass of the 
spherical dust configuration. The density p and the fluid 
velocity u diverge at any radius when t — > tc- Hence, we 
can simply solve for an expansion close to the center 



p(x) = p (t) + Pl (t)r + p 2 (t)r 2 + Oir 3 ) , 
u(x) = u it) + mit)r + u 2 it)r 2 + 0(r 3 ) . 



(6) 
(7) 



and the collapse occurs if these fields diverge at some t. 
The field equations impose pi(t) = Uo(t) = 1*2 (t) = and 

5 /3 

P2\t) = VPo (*)> wnere V < is a constant, and 



t(f>o) - t(pi) 



dxx 



V24^fG 



KT) 



3 1 ' / \ 

Pi + ^G {x - pih 



which, for k = 0, reduces to that in Ref. 0. This equa- 
tion can be integrated analytically for any n. A collapse 
occurs when tipo — > 00) > f(/9j), i.e. the time corre- 
sponding to an infinite density is in the future. We find 
that this condition is fulfilled only when k < but, when 
k > 0, the collapse does not occur. Fig. Q] shows that in 
this case the matter fields have an oscillatory behavior, 
whose period and amplitude depend on nr\. The same 
feature is found in early cosmology [l| . This suggests that 
singularities may be avoided in BF theory with k > 0, 
due to "repulsive gravity" effects proportional to up 1 . 
Ill A. Relativistic compact stars. Let us now con- 
sider static and spherically symmetric perfect fluid stars 
in the fully relativistic theory, described by 

q ab dx a dx b = ~pir)dt 2 + /i(r)dr 2 + r 2 dQ 2 , 
g ab dx a dx b = -Fir)dt 2 + S(r)rfr 2 + Air)r 2 d£l 2 . 

Here q a b is an auxiliary metric [l[ , and we have used the 
gauge freedom to fix the function in front of the spherical 
part of the metric q. We consider perfect-fluid stars with 
energy density p(r) and pressure P(r) such that 



T ab ^T^ r{ectfiuid = [p + P}u a u b + 



9 ab P, 



(8) 



where the fluid four- velocity u a = (f/vF,0,0,0) 
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We integrate the field equations (5) and (6) in Ref. [l[ 
imposing regularity conditions at the center of the star. 
The series expansion of the field equa tions at the center o f 
the star contains terms of the form — kP c )(1 + np c ). 
Assuming p Cl P c > 0, k must satisfy two conditions in 
order to allow for self-gravitating objects: 

P c k < 1 , for k > , (9) 
Pc\k\<1, for/t<0. (10) 

Hence, the existence of NSs with p c ~ 8 • 10 17 kg m -3 
and P c ~ 10 34 N m -2 strongly constrains the theory, 
\k\ < 1 m 5 kg _1 s -2 . Furthermore, it is easy to prove that 
compact objects only exist if P"(0) < 0. This gives a fur- 
ther constraint depending on p c , P c and p' c , whose form 
is cumbersome, but it is similar to Eq. ([5]). In particular, 
the condition is always satisfied for n > 0. 

The field equations are integrated outward up to the 
radius R, defined by the condition P(R) = 0, where 
we require the numerical solution to match the exact, 
and unique, vacuum Schwarzschild solution, F(r) = 
S(r) _1 = p(r) = ft(r) _1 = 1 - 2M/r, where M is 
the mass of the star. To match our numerically gen- 
erated spacetime to a Schwarzschild exterior we use the 
Darmois-Israel equations Q at the radius, i.e. [gij\ = 
and [Kij(q)] = 0, where [...] is the jump across the sur- 
face, Kij(q) is the extrinsic curvature tensor built with 
the metric q, and i,j = 0,2,3. These matching condi- 
tions come from the field equations and the requirement 
of a well-defined 3— geometry and give a unique prescrip- 
tion to compute the mass of the spacetime. 
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FIG. 2. Compactness for pressureless stars in the relativistic 
theory and in the Newtonian limit, as functions of np c . 

IIIA1. Relativistic pressureless stars. The existence of 
Newtonian pressureless stars makes it relevant to inves- 
tigate the existence of similar solutions in the full theory. 
To this purpose, we set P = 0. The conservation of the 
stress-energy tensor simply implies F(r) =const. The 
solutions of the field equations then depend only on one 
parameter, the dimensionless central density np c . 

As shown in Fig. [3J for any value of n > 0, there exists 
a regular solution which reduces to the Newtonian solu- 
tion discussed above in the nonrclativistic limit np c <C 1. 
These solutions have a positive binding energy and can 
be as compact as GM/R ~ 0.3 for np c ~ 200. Of course 



they do not exist in GR, while they exist in BF theory be- 
cause k > introduces a repulsive gravity contribution. 
Interestingly, the EOS for dark matter particles is ap- 
proximately P = 0. Hence, in this theory self-gravitating 
objects, purely made by dark matter, can exist and may 
reach the typical compactness of most compact NSs. Fur- 
thermore, these objects are stable in the Newtonian limit 
and it is reasonable to assume that they would remain 
stable also in the relativistic theory. 
Ill A 2. Polytropic EOS. We consider the model 

n m b ( n\ T ( n\ T 

p = nm b + K— — , P = Kn m b — 

F - 1 \n J V n o/ 

with the same polytropic parameters as in Ref. Q . Some 
results are shown in Fig. [3] for different values of k. The 
stellar mass M is shown as a function of the central bary- 
onic density p b — m b n(Q). In GR, maxima of this curve 
correspond to marginally stable equilibrium configura- 
tions, all solutions after the first maximum are unstable 
to radial perturbations (see e.g. (jf). This picture may 
change when k ^ 0. However, when np c -C 1 our so- 
lutions reduce to the non-relativistic ones, for which we 
proved stability, at least when k > 0. Hence, we conjec- 
ture that properties similar to GR still hold and branches 
before the first maximum in Fig. [3] arc likely stable. We 
leave a detailed analysis for the future. 

In the inset of Fig. [3] we also show the normalized 
binding energy reads E b /M = m/M — 1, where m = 
m b J d 3 Xy/—gu°n(r), is the baryonic mass of the config- 
uration and corresponds to the energy that the system 
would have if all baryons were dispersed to infinity. For 
bound (not necessarily stable) configurations, E b > 0. 

Positive values of k tend to enhance the relativistic 
effects: the maximum mass is larger than in GR and it 
occurs for smaller central density. Moreover, the binding 
energy for these models increases with k. Negative values 
of k have the opposite behavior. Remarkably, the most 
interesting effects show up when n > 0, i.e. in the same 
region where singularities seem to be prevented. 

These effects could be observable. Present NS obser- 
vations constrain the mass-radius relation (e.g. fioj). and 
electromagnetic observations of binaries containing X-ray 
pulsars may in principle constrain the binding energy as 
well [ll|. The recent discovery of a high-mass NS [l2| 
also rules out many EOS in GR. However, these observa- 
tions could be interpreted in terms of modified gravity at 
large curvature, rather then invoking exotic EOS in GR. 
IIIB. Slowly rotating models. Slowly rotating stars 
can be constructed from the corresponding static so- 
lutions (l3j . At first order in the rotation, g tip = 
—Q{r)r 2 sin 2 9, q tv = —r](r)r 2 sin 2 9 and the stress-energy 
tensor for a rotating fluid can be built from from Eq. © 
with 

«° = {« t J o,o J nu t } , «* = yf-(gtt + 2n tv + n 2 9 VV ), 

where f2 is the angular velocity of the fluid. The field 
equations for 77 and £ have to be solved by imposing reg- 
ularity at the center and matching the vacuum solution, 



4 




FIG. 3. Polytropic models for different values of k. Left panel: mass as a function of the central baryonic density p^. Right 
panel: mass-radius relation. Inset: binding energy as a function of pt- Results are normalized by po = 8 • 10 17 kg m~ 3 , which 
is a typical central density for NSs. Curves terminate when conditions Q or (110[) are not fulfilled. 



i] = C = 2J/r 3 at the stellar radius, where J is the angu- 
lar momentum. In Fig. |4]we show the moment of inertia 
I = J '/n as a function of the stellar mass. 
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FIG. 4. Moment of inertia for polytropic models as a function 
of the stellar mass for different values of «. 

IV. Conclusions. Eddington inspired theories are vi- 
able, one-parameter, alternatives to Einstein's gravity. 
We have shown that in these theories the structure of 
compact stars is dramatically different from GR, with 
potentially observable effects. For k > 0, our results 
show that BF theory has several remarkable features, 



e.g. singularities in gravitational collapse may be pre- 
vented. The mere existence of compact NSs strongly 
constrains the theory, kP c < 1. Furthermore, in our 
simple polytropic model, observational determination of 
the moment of inertia to an accuracy of 10%, as it is ex- 
pected from future observations of the double pulsar [lj| , 
will place even a stronger constraint, \npo\ < 0.1 (cf. 
Fig. 01). We expect that realistic EOS would constrain 
Kpo by the same order of magnitude. Remarkably, NSs 
are the densest matter configurations in the universe, so 
that these are likely the strongest bounds on the the- 
ory. Furthermore, it happens that the typical density 
of a NS, po ~ 8 • 10 17 kg m~ 3 , corresponds to the den- 
sity of the early universe (age ~ 10 _6 s), thus the present 
analysis can put strong constraints on the cosmological 
effects found in Ref. [lj. Several interesting issues, e.g. 
the collapse in the relativistic theory and black hole for- 
mation, the role of realistic EOS, the stability analysis 
of relativistic stars and possible crgoregion instability of 
rotating models, are left for future work. 
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